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Abstract 

We report evidence for the ground state of bottomonium, 775(16'), in the radiative decay T(3S) — )■ 
7% in e^e^ annihilation data taken with the CLEO III detector. Using 6 milhon T(35') decays, 
and assuming r(%) = 10 MeV/c^ we obtain B{T{3S) -/ijb) = (7.1 ± 1.8 ± 1.3) x 10"^, where 
the first error is statistical and the second is systematic. The statistical significance is ~4cr. The 
mass is determined to be M{r]b) = 9391.8 it 6.6 it 2.0 MeV/c^, which corresponds to the hyperfine 
splitting AMhfilS)b = 68.5 ± 6.6 ± 2.0 MeV/c^. Using 9 million T{2S) decays, we place an upper 
limit on the corresponding T(2S') decay, B{T{2S) jr]i,) < 8.4 x 10"** at 90% confidence level. 
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The spectroscopy of the bb bottomonium states provides valuable insight into Quantum 
Chromodynamics (QCD) since relativistic and higher-order corrections are less important 
for bb than any other qq system. Experimental measurements of the spectroscopic properties 
of the bottomonium states can therefore be compared with greater confidence with the 
predictions of perturbative QCD, as well as with lattice calculations. The hyperfine mass 
splitting of the singlet-triplet states is of particular interest since it probes the spin- dependent 
properties of the qq system. 

The triplet S state (I'^S'i) of bb bottomonium, T(15'), was discovered thirty years ago, 
but the identification of its partner, the singlet 5* state (1^5*0), rib{lS) (henceforth r/b), 
has eluded numerous searches, including those by CUSB [1], ALEPH |2J, DELPHI [3], 
and CLEO P]. As a result, the 15 hyperfine splitting, which is well-determined in the 
charmonium system, remained unknown in the bottomonium system. Recently, using 
their data sample of 109 million T(35') events, the BaBar collaboration reported E] 
the observation of the rji, with a statistical significance of more than lOa (standard de- 
viations) in the inclusive photon spectrum of T(3S') with the observed photon energy 
E^{T{3S) 'yr]b) = 921. 2t|^ ± 2.4 MeV, where the first error is statistical and the second 
is systematic. This gave M(r]h) = 9388. 9t2i ^ 2.7 MeV/c^ and a bottomonium hyperfine 
splitting, AMhf{lS)b = M(T(15)) - M{r]b) = 71.4l|^ ± 2.7 MeV/c^. BaBar's measured 
branching fraction was B{T{3S) — > 7?7b) = (4.8 ± 0.5 ± 0.6) x 10~^. Corroboration of the 
BaBar finding with an independent data set is essential. 

In this article we reexamine the CLEO data for the radiative decays T(3S', 2S) — > 777;,. An 
earher analysis of the same data resulted in upper hmits of B(T{3S) — )■ 777^) < 4.3 x 10~^ and 
B{T{2S) — )■ •yrib) < 5.1 x 10""^ at 90% confidence level [1]. However, the analysis had short- 
comings which are rectified in this article. The presence of the photon line corresponding to 
initial state radiation (ISR), e'^e~ — )■ 7T(1S'), located between the XbjC^P, 1-P) — ^ I'^i^S) 
region and the rjb signal region, was not included in the fits to the inclusive photon spectrum, 
an omission which biased the result toward small branching fractions. The assumption of 
^iVb) = MeV had a similar effect. Moreover, the analysis did not employ an important 
background-suppression variable, the angle between the radiative photon and the thrust axis 
of the rest of the event, introduced by BaBar [3]. We improve upon the previous publication 
by exploiting a more complete understanding of the expected photon line shape over a broad 
energy range to more accurately represent the Xfcj(2P, IP) — )• 7T(1S'), ISR, and rjb (with 
non-zero width) signals in a fit. We also employ a broader range of binning, fit ranges, and 
background parameterizations in order to avoid bias in any of these choices. 

The CLEO HI detector, which has been described elsewhere [7], contains a Csl elec- 
tromagnetic calorimeter, an inner silicon vertex detector, a central drift chamber, and a 
ring-imaging Cherenkov (RICH) detector, inside a superconducting solenoid magnet with a 
1.5 T magnetic field. The detector has a total acceptance of 93% of An. The photon energy 
resolution in the central (83% of Att) part of the calorimeter is about 2% at = 1 GeV 
and about 5% at 100 MeV. The charged particle momentum resolution is about 0.6% at 
1 GeV/c. 

The CLEO datasets correspond to (5.88±0.12) x 10*^ T{3S) and (9.32±0.19) x 10^ T{2S) 
decays. Our event selection for the inclusive photon spectra is identical to that reported 
in Ref. Events are required to have one or more photons, and three or more charged 
tracks. Photons with > 20 MeV are accepted in the "good barrel" region of the calorime- 
ter with |cos^| < 0.81 (where 6 is the polar angle with respect to the incoming positron 
direction), and are required to have a transverse spread which is consistent with that of 
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an electromagnetic shower. Photons from vr*^ decays are suppressed by vetoing any photon 
candidates that, when paired with another photon candidate in the good barrel or "good 
endcap" (0.85 < | cos6'| < 0.93) regions, have a mass within 2.5a of the known 7r° mass and 
cos 6*^^ > 0.7, where 9^^ is the opening angle of the photon candidates in the lab frame. 

We first consider the analysis of the inclusive photon spectrum from T{3S) decays. The 
analysis of T{2S) decays follows a similar path. In the region 500 < E^y < 1200 MeV, the 
spectrum consists of a peak centered around ^ 770 MeV due to the three unresolved 
transitions, XbjC^P) — ^ 7'^(1'S'), J = 0, 1, 2 on top of a smooth background that falls sharply 
with energy. The peaks due to ISR and rjb, which are more than an order of magnitude 
weaker than those from XbjC^P), are expected in the high energy tail region of the X6j(2P) 
peak. Hence, sensitivity to the possible presence of an rjfy signal depends critically upon 
properly representing the shape of the Xbji'^P) peaks as well as suppressing the underlying 
smooth background (as already achieved in part by the vr^ veto). As demonstrated by the 
BaBar analysis [5j, additional suppression can be achieved by recognizing that r]b signal 
photons are largely uncorrelated in direction with the rest of the event, whereas background 
photons from the continuum tend to follow the leading particles of the underlying event. 
This effect is more pronounced for T{3S) — )■ ■yrji, decays than for T(2S') — )■ jrib, but the effect 
is nevertheless useful for background suppression in both processes. The thrust angle {6t) is 
utilized to exploit these correlations; 6t is determined for each event as the angle between the 
momentum vector of the signal photon and the thrust vector [H] calculated using all other 
final state photons and charged particles boosted into the rest frame of the rji, candidate 
(defined by the signal photon). As shown in Fig. 1(a), the thrust angle distribution for the 
data events is peaked near | cos^tI = 1, whereas the thrust angle for the rjb signal events 
from Monte Carlo (MC) simulations is distributed uniformly. As a result, the sensitivity to 
a possible rjb signal in the presence of background varies greatly with | cos^tI, and it can be 
maximized by taking advantage of the | cos6't'| distribution. 

We utilize the |cos6'r| distribution, but in a manner quite different from that used by 
BaBar [5]. Instead of simply rejecting all events with large values of | cos 6t\, we increase the 
sensitivity to r/b by forming three separate photon energy spectra, one each for the | cos6't| 
regions (0.0,0.3) (I), (0.3,0.7) (II), and (0.7,1.0) (III), and performing a simultaneous joint 
fit to all three distributions. The signal-to-background ratio improves from region III to 
region II and from region II to region I, but all regions contribute to the sensitivity. Monte 
Carlo simulations show that, for a data sample of our size and a B{T{3S) — )■ 777b) whose 
value is assumed to be what is measured below, the three-region joint fit procedure leads to 
an average increase in the statistical significance ^ of an r]b signal of 0.6a over only accepting 
events with | cos < 0.7, albeit with a large r.m.s. spread of 0.7cr among MC trials. An 
average gain in significance over using no information about the thrust axis is 1.7a with 
an r.m.s. spread of 1.6a. Most of the 0.6a increase in sensitivity from the joint fit comes 
from splitting the | cos6'7^| < 0.7 region into two bins, which exploits the smaller background 
relative to expected signal in the | cos^^tI < 0.3 bin compared to the 0.3 < | cos^^tI < 0.7 
bin. On the average, inclusion of the | cos^^rl > 0.7 region by itself improves the result by 
0.2a. 

The photon peaks have shapes which are parameterized by convolving a relativistic Breit- 



We compute the statistical significance of the fit using the conventional likelihood expression 
■\/2 \n{Csig/Co), where Csig is the likelihood of the fit with a signal and Cq is the likelihood of the 
fit with the signal constrained to zero. 
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Wigner resonance function with a Crystal Ball (CB) calorimeter response function [9], which 
consists of a Gaussian part with width a (the energy resolution) smoothly joined to a low- 
side power-law tail described by two additional shape parameters. The energy resolution and 
CB shape parameters were determined with two complementary methods. In Method A, we 
utilized isolated photons in e'^e~ — )■ e'*"e~7 data events with photon energies near -Etme = 
750 MeV, where -Etrue is the photon energy expected from using only the measured angles of 
the and 7. We then extracted an inherent line shape by deconvolving the spread in -Etrue 
(obtained from simulated events) from the observed E^/Etme- In Method B, we compared 
exclusive T(35) ^ 7Xw(2P), Xbi{2P) 7T(15), T{1S) £+r (£± = e± or /i±) in data 
and MC simulation to determine the shape of the T(35') — )■ 7X;,i(2P) photon line. The data 
distribution was used to determine the Gaussian part of the shape and the MC simulations 
were used to determine the two tail parameters after tuning the MC parameters to match the 
Gaussian part observed in the data. Methods A and B lead to consistent energy resolutions 
and CB shape parameters, resulting in a line shape that is significantly different from that 
used in the original CLEO analysis. While the tail parameters of the peak shapes are fixed 
to be the same for all three relevant photon energies (xbj(2P), ISR, and rjb), the Gaussian 
widths for the three are different. The fitted Gaussian width for the overlapping Xfej(2P) 
peaks near 770 MeV in the inclusive spectrum is ct(770 MeV) = 16.7 ± 1.0 MeV. The 
variation of the photon resolution width with energy was determined from MC simulations 
made for a wide range of photon energies. Its parametrization was used to obtain the 
extrapolated values, cr(859 MeV) = 17.4 ± 1.0 MeV, and a(920 MeV) = 18.3 ± 1.1 MeV, for 
the ISR and 77^ peaks, respectively. 

The expected intensity of the ISR peak was obtained by extrapolating its yield observed 
in CLEO data taken on the T(4S') resonance. The expected yield iV(ISR) = 1726 ± 131, 
photon energy ii^^(ISR) = 859 MeV, and energy resolution (t(ISR) = 17.4 MeV are fixed in 
all fits of the inclusive spectra. 

The prominent peaks in the inclusive spectra shown in Fig. 1(b) are composites of the 
three X6j(2P) — ?■ 7T(15') peaks for J = 0,1,2. We fix the relative strengths of these 
three lines to the ratios determined from other measurements [10] and float only the overall 
amplitude. We also fix the spin-orbit splitting of these lines to the values measured in 
Ref. [1], but we float the absolute energy scale. The latter provides a useful check on our 
uncertainty in the absolute energy calibration. The CB line shape parameters are fixed as 
discussed previously, while the effective energy resolution, which includes Doppler smearing, 
is allowed to float. 

The efficiencies for Xbj(2P), ISR, and rji, in our event selections are obtained by Monte 
Carlo simulations with the 1 + a cos^ 6 angular distributions expected for El and Ml tran- 
sitions with appropriate values of a for X6i(2P) and Xfe2(2P), and a = 1 for XboC^P) and 
rif,. Separate calculations were done for the three | cos^^l bins, and it was found that effi- 
ciencies are approximately constant in | cos^t|. The summed efficiencies for r]), and ISR are 
(54.2 ± 3.8)% and (6.9 ± 0.1)%, respectively. 

As discussed previously, we perform a joint fit of the data in three | cos 6t\ bins. All fitting 
parameters (apart from those in the background function described below) are constrained 
to be the same in the three | cos6't| bins. That is, the yields for the Xftj(2P), ISR, and rjb 
photon peaks in each of the three | cos^tI bins were constrained to be proportional to the 
ratios A| cosOxli/^i where is the signal efficiency for bin i. 
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The smooth background was fitted with exponential polynomials, 

j^^exp(g<,E;). (1) 

As the only experimental handle on these backgrounds is the inclusive spectrum itself, we 
explored uncertainties in their determination by varying binning types (both linear and 
logarithmic binning were used), the order of the polynomial {n was varied from 2 to 4 in 
each thrust bin independently) and the fit range (six different ranges were tried extending 
down to 500 MeV and up to 1340 MeV). Results for the rji, (mass, significance, and branching 
fraction) were then averaged through all fits with confidence level (CL) above 10%. The 
r.m.s. spread among the fit variations was taken as a measure of the systematic uncertainty in 
the background determination. Averaged through all successful fits, the maximum likelihood 
significance of the rji, signal is 4. la with a r.m.s. spread of 0.4. A representative fit, whose 
parameters are close to the average values for the ensemble of accepted fits, is chosen as 
nominal. This fit (shown in Fig. 2) has N{'yr]i,) — 2311 ± 546 counts and gives B{'y7]b) = 
B{T{3S) 7776) = (7.1 ± 1.8) x 10"^ and E^{^r)b) = E^{T{3S) -frjb) = 918.6 ± 6.0 MeV, 
with a CL of 18.5%. 

TABLE I: Summary of estimated systematic uncertainties and their sums in quadrature for the 
T{3S) 7?7(, analysis. The item labeled Background refers to variations of the background function 
parameters, the fit range, and linear versus logarithmic binning. 



Uncertainty in 



Source 


£;^(7%) (MeV) 


B{im) (%) 


Background 


±1.0 


±10 


Photon Energy Calibration 


±1.2 




Photon Energy Resolution 


±0.3 


±2 


CB and Xhj{'^P) Parameters 


±0.7 


±8 


ISR Yield 


±0.4 


±3 


Photon Reconstruction 




±2 


iV(T(35)) 




±2 


MC Efficiency 




±7 


i]h Width 


±0.6 


±9 


Quadrature sums 


±1.9 


±18 



The systematic uncertainties in our results are obtained as follows and are summarized 
in Table 1. We assign the r.m.s. variations in the results obtained for all the accepted fits, 
±1.0 MeV in E^{'~frii,), and ±10% in Bi^-^rjiy) as systematic uncertainties due to background 
shape, binning, and range variations. The changes in our results are negligible when we 
alter the lower CL hmit for acceptable fits from 10% to either 5% or 15%. We vary the 
photon energy resolution, the Crystal Ball shape parameters, and the Xbji'^P) parameters 
within their errors and assign the resulting variations in E^{'^r]b) and B{'yr]b) as systematic 
uncertainties. 

The fitted Xbj(2P) centroid energy in our data is 769.9 ± 0.2 MeV, while the expected 
energy is 769.6l?:S MeV. The 0.3 MeV deviation of our measured value suggests that 
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our photon energy calibration has a maximum possible uncertainty of 1?;2 MeV. This is 
consistent with our measurement of ISR photon energies from T(45') and below T(4S') data, 
which agree with the expected energies within ±0.3 MeV. Based on these considerations 
we conservatively assign the systematic uncertainty due to photon energy calibration as 
±1.2 MeV. We obtained the value of Bi^rjb) by assuming r(?7;,) = 10 MeV/c^. We find 
that B{-yrih) depends linearly on the assumed value of T{rih) in MeV/c^, as Bi^jrib) = [5.8 + 
0.13 r(?7;,)] X 10~^. Varying r(?7b) from 5 to 15 MeV/c^, a range that includes nearly all 
theoretical expectations [11], the branching fraction changes by ±0.65 x 10""^ or ±9%. This 
uncertainty in the rji, width also contributes ±0.6 MeV to the uncertainty in E^{'~frii,). Other 
systematic uncertainties are due to the Monte Carlo efficiency calculation and the number 
of T(3S') events. 

In fitting the 777;, peaks, we do not include the factor [12] {E^/Eq)'^[1 + a{E^/Eo)'^]'^ 
expected in the decay width for the hindered Ml transition T(35') — >■ jrib. {Eq is the photon 
energy for the central value of the rjb mass.) While theoretical estimates [12] of alpha vary, 
a = 1 leads to a distortion of the rjh peak shape and a consequent reduction of E^ipfTjb) by 
approximately 3 MeV. Since our data sample is not large enough to determine a, in the 
absence of firm theoretical predictions we do not include this effect as a bias or as a term in 
our systematic error. 

Our final results are: ^^(7//;,) = 918.6±6.0±1.9 MeV and (77^6) = (7.1±1.8±1.3) x 10-^ 
where the first errors are statistical and the second errors are systematic. Our result for 
E^(7?7b) corresponds to M(r/f,) = 9391.8 ± 6.6 ± 2.0 MeV/c^ and ^Mhf{lS)b = 68.5 ± 6.6 ± 
2.0 MeV/c^. This is consistent with lattice QCD predictions that employ dynamical quarks 
and include both continuum and chiral extrapolations [13]. Our results for both AMhf{lS)j, 
and B{'~frih) are also well within the wide range of pQCD based theoretical predictions [13] ■ 
Both measurements are in good agreement with the BaBar measurements |[5l |6j . 

We also analyzed our data set containing (9.32±0.19) x 10^ T(2S') events for T(2S') — )■ •yrib 
using the same event selection and joint fit analysis procedure as described above for 
T(3S') — )■ '-yrib. One difference is that we chose to represent the T(2S') — )■ 7r°7r°T(15') back- 
ground component explicitly in the fit since it introduces a kink in the spectrum not far from 
the signal region. The shape of this background was taken from Monte Carlo simulations. 
Its normalization was fixed to the PDG value of the branching fraction. Unlike in the T(25') 
analysis, the addition of the explicit T(3S') — tt^tt^T{1S) background component to the 
T(3S') fits had a negligible effect on the results. In the expected signal region for T(2S') ra- 
diative decay, 200 < E^ < 800 MeV, the background is an order of magnitude larger than in 
the T(3S') signal region, and in none of the T(2S') | cos 9t\ regions could the r]h be identified. 
In the joint fit analysis (shown in Fig. 3), fixing ii^^(T(2S') — )■ 777;,) = 611 MeV, correspond- 
ing to rjf, mass determined in T(35') decay, leads to ;B(T(25') — )■ 777b) = (3.9±2.7±2.3) x 10~^, 
or an upper limit of ;B(T(25') — )■ 777b) < 8.4 x 10~^ at 90% confidence level. This is consistent 
with the BaBar T(25) result [6], B{T{2S) 77/;,)) = {3.9til ± 0.9) x 10"^ 
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FIG. 1: (a) Distribution of | cos 0^1 for r]{y signal MC events (dotted) and background dominated 
T{3S) data (shaded) in three regions, I, II, and III defined in the text. The histogram of the 
Monte Carlo simulation of the T(3S') — )• 777^ signal has been multiplied by a factor of 3000 to make 
it visible, (b) The distribution from T(35') data in the three regions of |cos^t|- Only the 
XbjC^P) — ^ 7^(15) lines at around 770 MeV are visible above the background. 
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FIG. 2: Background and X6j(2-P) subtracted distributions of E-./ from T{3S) decays in three | cos 9t\ 
regions, I, II, and III defined in the text. The curves are the results of the joint fit, with a CL of 
18.5%. The Xfej(2-P) peaks are indicated by the dotted Hues and the rjb signals by the dashed lines, 
which join the solid line. 
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FIG. 3: Background and Xbji^P) subtracted distributions of from T{2S) decays in three | cos 9t\ 
regions, I, II, and III defined in tlie text. The curves are the joint fit results. The Xbj{^P) peaks 
are indicated by the dotted hues and the 90% rjb upper hmits by the dashed hues. 
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